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Abstract. The ground-state magnetic phase diagram of a spin S — 1/2 two-leg ladder with alternating 
rung exchange J_L(n) = Jj_[l -|- {—1)"S] is studied using the analytical and numerical approaches. In the 
limit where the rung exchange is dominant, we have mapped the model onto the effective quantum sine- 
Gordon model with topological term and identified two quantum phase transitions at magnetization equal 
to the half of saturation value from a gapped to the gapless regime. These quantum transitions belong to 
the universality class of the commensurate-incommensurate phase transition. We have also shown that the 
magnetization curve of the system exhibits a plateau at magnetization equal to the half of the saturation 
value. We also present a detailed numerical analysis of the low energy excitation spectrum and the ground 
state magnetic phase diagram of the ladder with rung-exchange alternation using Lanczos method of 
numerical diagonalizations for ladders with number of sites up to A'^ = 28. We have calculated numerically 
the magnetic field dependence of the low-energy excitation spectrum, magnetization and the on-rung spin- 
spin correlation function. We have also calculated the width of the magnetization plateau and show that 
it scales as 5", where critical exponent varies from v — 0.87 ± 0.01 in the case of a ladder with isotropic 
antiferromagnetic legs to v = 1.82 ±0.01 in the case of ladder with ferromagnetic legs. Obtained numerical 
results are in an complete agreement with estimations made within the continuum-limit approach. 

PACS. 75.10.Jm Quantized spin models; 75.10.PqSpin chain models 



1 Introduction 



Low-dimensional quantum magnetism has been the sub- 
ject of intense research activity since the pioneering pa- 
per by Bethe [1] . Perpetual during the decades interest in 
study of these systems is determined by their remarkably 
rich and unconventional low-energy properties (see for re- 
view Ref. 2). An increased current activity in this field is 
connected with the large number of qualitatively new and 
dominated by the quantum effects phenomena discovered 
in these systems OH] as well as with the opened wide 
perspectives for use low-dimensional magnetic materials 
in modern nanoscale technologies. 

The spin S = 1/2 two- leg ladders represent one, partic- 
ular subclass of low-dimensional quantum magnets which 
also has attracted a lot of interest for a number of rea- 
sons. On the one hand, there was remarkable progress in 
recent years in the fabrication of such ladder compounds 
[5]. On the other hand, spin-ladder models pose interest- 
ing theoretical problems since antiferromagnetic two-leg 
ladder systems with spin 5=1/2 have a gap in the ex- 
citation spectrum [6ll7l[8] they reveal an extremely rich 
behavior, dominated by quantum effects in the presence 
of a magnetic field. These quantum phase transitions were 



intensively investigated both theoretically [9-24] and ex- 
perimentally [25-30]. 

Usually, these most exciting properties of low dimen- 
sional quantum spin systems exhibit strongly correlated 
effects driving them toward regimes with no classical ana- 
log. Properties of the systems in these regimes or "quan- 
tum phases" depend in turn on the properties of their 
ground state and low-lying energy excitations. Therefore 
search for the gapped phases emerging from different sources 
and study of ordered phases and quantum phase tran- 
sitions associated with the dynamical generation of new 
gaps is an important direction in theoretical studies of 
quantum spin systems. 

A particular realization of such scenario appears in 
the case where the spin-exchange coupling constants are 
spatially modulated. The spin-Peierls effect in spin chains 
represent prototype example of such behavior [31| . In the 
recent paper of one of the authors the new type of effec- 
tive spin-Peierls phenomenon in ladder systems, connected 
with spontaneous dimerization of the system during the 
magnetization process via alternation of rung exchange 
has been discussed [35] • The Hamiltonian of the model is 
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+ Ji^[l + (-l)"(5]S„a •S„,2, (1) 

n 

where Sn.a is the spin S = 1/2 operator of rung n (n=l,...,L) 
and leg a (a = 1, 2). The interleg coupHng is antiferromag- 
nctic, Jf = J±{1±5) > 0. 

In Rcf. 32 the model was studied analytically in the 
limit of strong rung exchange and magnetic field ~ 
H ^ |J|||,(5Jj^. In this limit the ladder Hamiltonian is 
mapped onto the spin- 1/2 XX Z Heiscnbcrg chain in the 
presence of both longitudinal uniform and staggered mag- 
netic fields, with the amplitude of the staggered compo- 
nent of the magnetic field proportional to ~ 6J^. The 
continuum-limit bosonization analysis of the effective spin- 
chain Hamiltonian show, that the alternation of the rung- 
exchange leads to the dynamical generation of a new en- 
ergy scale in the system and to the appearance of two addi- 
tional quantum phase transitions in the magnetic ground 
state phase diagram. These transitions manifest them- 
selves most clearly in the presence of a new magnetization 
plateau at magnetization equal to one half of its satura- 
tion value. It was shown that the new commensurability 
gap (correspondingly the width of magnetization plateau) 
scales as S^, where f = 4/5 in the case of a ladder with 
isotropic antiferromagnetic legs and ^ = 2 in the case of 
a ladder with isotropic ferromagnetic legs. Although up 
to now no materials are available which realize this mod- 
els, theoretical studies of this type of Peierls distortion is 
extremely interesting, since such an instability could in 
principle develop during the magnetization process of an 
antiferromagnetic ladder, in particular with high proba- 
bility in the case of applied uniaxial along legs pressure. 

In this paper we continue our studies of an isotropic 
spin 5* = 1/2 two-leg ladder with alternating rung ex- 
change in a uniform magnetic field using the numerical 
analysis. In particular we apply the Lanczos method to 
diagonalize numerically finite ladder systems with lengths 
L = 6,8,10,12,14. Using the exact diagonalization re- 
sults, we calculate the spin gap, magnetization and the 
intra-rung spin correlations as a function of applied mag- 
netic field. We have also calculated the spin-density dis- 
tribution in the ground state at magnetization plateau. 
Based on the exact diagonalization results we obtain the 
ground-state magnetic phase diagram of the model show- 
ing four quantum phase transitions, in agreement with the 
predictions made in Ref. 32. We also numerically com- 
puted the width of the plateau and show it scales as S'^ , 
where v = 0.87±0.01 in the case of a ladder with isotropic 
antiferromagnetic legs and i/ = 1.82 ± 0.02 in the case of 
a ladder with isotropic ferromagnetic legs. 

The paper is organized as follows. In the forthcom- 
ing section we briefly summarize results of the analytical 
studies. In section III, we present numerical results of our 
exact diagonalization studies of the system. Finally, we 
conclude and summarize our results in section IV. 

2 Effective Hamiltonian 

In this section we briefly summarize the results obtained 
within the analytical studies. To obtain the effective spin- 



chain we follow the route already used to studies of the 
standard two-leg ladder models in the same limit of strong 
rung exchange [T21IT3] . We start from the case J|| = 0, 
where the system decouples into a set of noninteracting 
rungs with couplings J_\_{n) — J±\l ± (— 1)"(5] and an 
eigenstate of the Hamiltonian is written as a product of 
rung states. At each rung, two spins form cither a singlet 
state \s) with energy Es{n) = — 0.75Jj_(ri) or in one of 
the triplet states \t'^),\t^) and \t^) with energies Ei+{n) = 
0.25J±{n)-H,Eto{n) = 0.25 J±{n) and Et- (n) = Q.25J±{n) + 
H, respectively. When H is small, the ground state con- 
sists of a product of rung singlets. As the field H in- 
creases, the energy of the triplet state 1^+) decreases and 
at _ff ~ J±{'n) forms, together with the singlet state, a 
doublet of almost degenerate low energy state, split from 
the remaining high energy two triplet states. This allows 
to introduce the effective spin operator r which act on 
these states as [12] 

^,^|S0>„ = \t+>n, T+\t+ >„ =0, (2) 
T~\so>n = 0, T'lf^ >n= \so >n ■ 

The relation between the real spin operator S„ and the 
pseudo-spin operator t„ in this restricted subspace can 
be easily derived by inspection. 




Using Q, to the first order and up to a constant, we easily 
obtain the effective Hamiltonian 

^e// "^{^J^y {Tn'Tn+i + h.C.) + Jz«+i} 
n 

-/^e//E<-^e//E(-l)"<' (4) 

n n 

where Jxy = 2J^ = Jy and 

h1^f^H-J^-J\^/2, hl^f^SJ^. (5) 

Thus the effective Hamiltonian is nothing but the XXZ 
Heisenberg chain, with anisotropy Jz/ Jxy = ^ = 1/2 in a 
uniform and staggered longitudinal magnetic fields. It is 
worth to notice that a similar problem has been studied 
intensively in recent years [551IM1I55 11361157] . 



2.1 The first critical field i?ci and the saturation field 



The performed mapping allows to determine critical fields 
TJci corresponding to the onset of magnetization in the 
system and the saturation field [12]. The easiest way 
to express Hc^ and in terms of ladder parameters is 
to perform the Jordan- Wigner transformation which maps 
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the problem onto a system of interacting spinlcss fermions 
[38]: 

Hsf = i ^(a+a„+i + h.c.) + V'^pnPn+i 

n n 



where t ~ V ~ J\\l'^i Mo = h^eff + "^11 and /^i = h\^j. 
The lowest critical field Hc^ corresponds to that value of 
the chemical potential /Ltoc for which the band of spin- 
less fermions starts to fill up. In this limit we can ne- 
glect the interaction term in Eq. ^ and obtain the model 
of free massive particles with spectrum E^{k) = — /io ± 

(J2c0s2(fc) +^2)l/2_ rpj^ig gives = J_l-(J|+5Vi)l/2. 

A similar argument can be used to determine ■ In 
the limit of almost saturated magnetization, it is useful 
to make a particle-hole transformation and estimate 
from the condition where the transformed hole band starts 
to fill, what gives i?c2 = -^J- 



(j2+52ji)l/2. 



2.2 Magnetization plateau: 

To determine parameters characterizing the magnetiza- 
tion plateaux at M = O.SMgat, we use the continuum- 
limit bosonization treatment of the model ([4]). To obtain 
the continuum version of the Hamiltonian (j4]) we use the 
standard bosonization expression of the spin operators [39] 



f— 1)" 

d,:(t) + — ain{VATTK(j)) , (7) 




2tt 



(-1)" + sin\/47rA'<?!) 



(8) 



where (l){x) and 0{x) are dual bosonic fields, dt4' = v^OxO 
and the spin-stiffness parameter K = Tr/2 (tt — arccosZi). 
Using ((Zj)-® we get the following bosonized Hamiltonian 

1(^4^0) I (9) 



TTflo 



sm 



which is easily recognized as the standard Hamiltonian for 
the commensurate-incommensurate phase transition [40] 
and the Bethe ansatz ^l] ■ We use these results to describe 
the magnetization plateau and transitions from a gapped 
(plateau) to the gapless Luttinger liquid phase. 

Let us start our consideration from the case h'^fj = 
corresponding to H = J± + J\\/2, where the Hamiltonian 
^ reduces to the quantum sine-Gordon (SG) model with 
a massive term ~ /i^yj sin(v47r^(/)). From the exact so- 
lution of the SG model [42] it is known that the excita- 
tion spectrum is gapfuU for0<i4r<2. Atl<isr<2 
the excitation spectrum of the model consists of solitons 
and antisolitons with mass M, while for < A' < 1 the 



spectrum contains also soliton-antisoliton bound states 
("breathers"). However for 1/2 < X < 1 the soliton mass 
M remains the lowest excitation energy scale in the model. 
Below we will restrict our consideration by limiting cases 
of the isotropic antiferromagnetic and ferromagnetic lad- 
der, corresponding to the case Z\ = 1/2 [K = 3/4) and 
A = —1/2 [K = 3/2), respectively. Thus, in both con- 
sidered cases the effective bosonized Hamiltonian is in the 
gapped regime with gap determined by the soliton mass 
M . Moreover from the exact solution [43] we can also get 
the exact scaling relation between the soliton mass phys- 
ical M and it's bare value toq = 5J± 



1/{2~K) 



(10) 



Thus at hP^ff = the excitation spectrum of the system 
is gapped. In the ground state the field 4> is pinned in one 
of the minima where (01 sin (V47rA»|0) = -1, what cor- 
responds (see Eq. (O) to a long-range antiferromagnetic 
order in the ground state of the effective Heisenberg chain, 
i.e. to a phase of the initial ladder system, where odd rungs 
have a dominant triplet character and even rungs are pre- 
dominantly singlets. 

At h^i.ff 7^ the very presence of the gradient term 
in the Hamiltonian ^ makes it necessary to consider the 
ground state of the sine-Gordon model in sectors with 
nonzero topological charge. The effective chemical poten- 
tial ~ h'^fj:dx4> tends to change the number of particles 
in the ground state, what immediately implies that the 
vacuum distribution of the field (f> will be shifted with re- 
spect of the minima which minimize the staggered part. 
This competition between contributions of the smooth 
and staggered components of the magnetic field is re- 
solved as a continuous phase transition from a gapped 
state at \h^ff\ < M to a gapless Luttinger liquid phase 

at |/ie//l > where M is the soliton mass [40]. This 
condition immediately gives two additional critical values 
of the magnetic field 



H± = J, 



,/||/2±M, 



(11) 



and respectively the width of the magnetization plateau 
given by 



H+ - h: 



2J|| (<5J_i/J|| 



>l/(2-if) 



(12) 



It is straightforward to get, that in the case of a ladder 
with antiferromagnetic legs, the width of the magnetiza- 
tion plateau scales as S^^^ while in the case of a ladder 
with ferromagnetic legs as S^. 

In order to investigate the detailed behavior of the 
ground state magnetic phase diagram and to test the va- 
lidity of the picture obtained from continuum-limit bosoniza- 
tion treatment, below in this paper we present results of 
numerical calculations using the Lanczos method of exact 
diagonalizations for finite ladders. 
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3 Numerical results 



In this section, to explore the nature of the spectrum and 
the quantum phase transition, we used Lanczos method to 
diagonalize numerically ladders with length up to L = 14. 



3.1 The Energy Gap 

First, we have computed the three lowest energy eigenval- 
ues of ladders with antiferromagnetic (ferromagnetic) legs 
J|| = 1.0(— 1.0) and different values of the rung exchanges 
J±{n). To get the energies of the few lowest eigenstates 
we consider ladders with periodic boundary conditions on 
legs. 

In FiglU we present results of these calculations for 
the rung exchanges J° = 11/2,6 = 1/11 and ladder sizes 
L = 8, 10. We define the excitation gap as a gap to the 
first excited state. As it can be seen in FigdJ in the con- 
sidered limit of strong on-rung exchange, this difference 
is characterized by the indistinguishable (within the used 
numerical accuracy) dependence on the ladder length and 
shows an universal linear decrease with increasing mag- 
netic field. 

At _ff = the spectrum model is gapped for all cases 
(ferromagnetic and antiferromagnetic legs) . For H ^ the 
energy gap decreases linearly with H and vanishes at Hc^ . 
This is the first level crossing between the ground-state 
energy and the first excited state energy. The spectrum 
remains gapless for < H < H~ and once again be- 
comes gapped for H > . The spin gap, which appears 
at H > , first increases with the external uniform mag- 
netic field, but then starts to decrease, and again vanishes 
at H^. With more increasing field H > the spec- 

trum remains gapless up to the critical saturation field 
Hc2- Finally, a,t H > Hc.^ the gap opens again and for a 
sufficiently large field becomes proportional to H. In the 
regions Hc^ < H < H~ and < H < Hc2, the two 
lowest states cross each other N/2 times, which causes 
gapless regimes in the thermodynamic limit N — > oo. 
These level crossings lead to incommensurate effects that 
manifest themselves in the oscillatory behavior of the re- 
sponse functions. 

Let us now comment the numerical method which we 
use to extract critical values of the magnetic field in the 
thermodynamic limit (L —f oo) from the finite ladder ex- 
act diagonalizations data. In an infinite system the critical 
point corresponds to a value of the magnetic field at which 
the energy gap vanishes. Since, in the case of finite systems 
the gap remains finite, we use the phenomenological renor- 
malization group (PRG) method (see for details Ref.44 ) 
to obtain values of the critical fields. The phenomenolog- 
ical renormalization group equation for the energy gap 
G{L, H) is written as 



[L + 2)G{L + 2, H') = LG{L, H) . 



(13) 



At the quantum critical point, LG{L, H) should be size 
independent for large enough systems in which the con- 
tribution from irrelevant operators is negligible. This al- 
lows to determine critical points very accurately. We define 
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Fig. 1. Difference between the energy of the first excited state 
and the ground state energy as a function of the magnetic field 
H, for ladder with rung exchange J± = 11/2, S = 1/11 and 
with (a) ferromagnetic legs J\\ — —1.0 (b) antiferromagnetic 
legs Jii — 1.0, including different ladder lengths L = 8, 10. 



Hc{L, L + 2) as L-dependent fixed point of PRG equation 
and then it is extrapolated to the thermodynamic limit in 
order to estimate He- At the first step we plot the curve 
LG{L, H) versus H for system sizes L and L + 2. These 
curves cross at a certain value Hc{L, L + 2) which is deter- 
mined as a finite size critical point. The thermodynamic 
critical point He is obtained by extrapolating Hc{L, L + 2) 
to L — > oo. Figure [2] shows the extrapolation procedure 
of the transition points for the ladder with antiferromag- 
netic legs and rung-exchange parameters J± = 11/2 and 
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Fig. 2. Fixed points Hc{L,L + 2) of the PRG equation 
for the energy gap, plotted versus 1/iL + 1), for Jy = 1.0, 
J_L = 11/2,5 = 1/11 and different ladder lengths L — 
4,6,8, 10, 12, 14. The values otHc{L,L + 2) for five pairs of sys- 
tem sizes (L,L + 2) = (4, 6), (6, 8), (8, 10), (10, 12), and (12,14) 
are represented in the figure. 



S = 1/11. Obtained by this method vahies of critical are 

= 4.48 ±0.01, i7c2 = 7.62 ±0.01, (14) 

(15) 



H- = 5.32 ±0.01, 



H+ = 6.87 ±0.01. 



Analogously in the case of the ladder with ferromagnetic 
legs we obtain the following set of critical parameters 



4.47 ±0.01, 



H- = 4.89 ±0.01, 



H. 



5.62 ±0.01, (16) 
5.26 ±0.01. (17) 



One can easily check that the obtained values for critical 
field obtained from the finite ladder studies are very close 
to the values predicted analytically. 



3.2 Magnetization curve 

To study the magnetic order of the ground state of the 
system, we start with the magnetization process. First, 
we have implemented the Lanczos algorithm on finite lad- 
ders to calculate the lowest cigenstate. The magnetization 
along the field axis is defined as 



1 ^ 

jY.^Gs\{Sl^ + Sl,,)\Gs), 



(18) 



n=l 



where the notation {Gs\...\Gs) represent the ground state 
expectation value. In Fig. [31 we have plotted NF- as a 
function of the magnetic field i7, and for a ladder with 
(a) ferromagnetic and (b)antiferromagnetic legs and with 



0.5 - 



0.4 - 



0.3 - 



ANTIFERROMAGNETIC 
LEGS 

MAGNETIZATION PLATEAUX 




0.2 - 



Fig. 3. a. The magnetization along field as a function 
of the applied magnetic field H for (a) ferromagnetic ladder 
J|l = —1.0 (b) antiferromagnetic ladder Jy = 1.0, including 
different ladder lengths L = 6, 8, 10, 12. 



rung exchange parameters = 11/2,(5 = 1/11 and dif- 
ferent lengths L = 6, 8, 10, 12 . 

As it is clearly seen in Fig. [3] besides the standard 
rung-singlet and saturation plateaus at H < Hc^ and 
H > respectively observed in the case of ladders with 
uniform rung exchange pillOj. we observe a plateau at 
M = 0.5AIsat- Observed oscillations of the magnetization 
at < H < H~ and iJ+ < H < Hc^ result from the 
level crossing between the ground and the first excited 
states of this model in the gapless phases. Usually, in or- 
der to give an estimation of the width of magnetization 
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ization order parameter) for even and add sites, as 



0.25 
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Fig. 4. The on-rung spin correlation functions for odd (J^ — 
5.0) and even ( — 6.0) rungs as a function of the applied 
field H for L = 8, 10 lengths with rung exchange parameters 
J± = 11/2,5 = 1/11 with (a) ferromagnetic (Jy — —1.0) and 
(b) antiferromagnetic legs (Jy — 1.0) . 



plateau in the thermodynamic limit, the size scaling of its 
width is performed [H] . We also perform the similar anal- 
ysis for selected values of the exchange parameters and 
found that the size of the plateau interpolates to finite 
value with L — > oo. 

An additional insight into the nature of different phases 
can be obtained studying the intra-rung correlations. We 
define the on-rung spin correlation function (rung dimer- 



dt^-J2(Gs\S,,2m-S2,2m\Gs) 



and 



2m+l 



^+l\Gs) 



(19) 



(20) 



taking sum over even or odd sites, respectively. In Fig. S] 
we have plotted the dj; and d° as a function of the mag- 
netic field H for the ladder of lengths L = 8, 10 with (a) 
ferromagnetic and (b) antiferromagnetic legs and the rung 
exchange parameters J± = 11/2,(5 = 1/11. As it is seen 
from this figure, at H < H^^ spins on all rungs are in a 
singlet state = d° ~ —0.75, while a,t H > Hc^, dr is 
slightly less than the saturation value = d° ~ 1/4 and 
the ferromagnetic long-range order along the field axis is 
present. Deviation from the saturation values —3/4 and 
1/4 reflects the effect of quantum fluctuations. However, 
in the considered case of strong rung-exchanges and high 
critical fields quantum fluctuations are substantially sup- 
pressed and calculated averages of on-rung spin correla- 
tions are very close to their nominal values. 

On the other hand, for intermediate values of the mag- 
netic field, at i?ci < H < Hc^ the data presented in Fig. 2] 
gives us a possibility to trace the mechanism of singlet- 
pair melting with increasing magnetic field. As it follows 
from Fig. |4] at slightly above iJci spin singlets pairs 
start to melt in all rungs simultaneously and almost with 
the same intensity. This effect is more profound in the case 
of ladder with ferromagnetic legs (Fig. Hja) where the on- 
rung spin correlation functions increase almost parallel for 
even and odd rungs. With further increase of H melting 
of weak rungs gets more intensive, however at H = 
the process of melting stops. As it is seen in Fig. IHb weak 
rungs are polarized, however their polarization is far from 
the saturation value d° ~ 0.1, while the strong rungs still 
manifest strong on-site singlet features with d'^ ~ —0.62. 
This effect is much weaker in the case of ferromagnetic 
legs in accordance with previous results. Moratorium on 
melting stops at H = , however for H > strong 
rungs start to melt more intensively while the polariza- 
tion of weak rungs increases slowly. Finally at H = Hc^ 
both subsystems of rungs achieve an identical, almost fully 
polarized state. Note, that the almost symmetric fluctua- 
tions in on-rung correlations, increase in d'^ at H < H~ 
decrease in d° at H > reflect the enhanced role of 
quantum fluctuations in the vicinity of quantum critical 
points. 

To complete our description of the phase at magneti- 
zation plateau with M = O.bMsat we have calculated the 
rung-spin distribution in the ground state 



(21) 



In Fig. O we have plotted the spin distribution in the 
ground state of a ladder with rung-exchange parameters 
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-0.4 - 



-0.6 - 



Fig. 5. The spin distribution in the GS as a function of the 
rung number "n" for magnetization corresponding to plateau 
at — 0.5Mgat and for rung exchanges Jj^ = 5.0 and — 
6.0. 
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-1.6 



-1.4 



-1.2 



log (5) 



J± = 11/2, (5 = 1/11 and antiferromagnetic legs as a func- 
tion of the rung number "n" for a value of the magnetic 
field corresponding to the plateau at = 0.5M|^j. To 
obtain an accurate estimate of the function M^, we have 
calculated them from the Eq.( [2T|) for system sizes of 
L = 6,8,10,12,14. The thermodynamic limit {L oo) 
of the finite size results are obtained by extrapolation 
method and used for plotting. As we observe the sys- 
tem shows a well pronounced modulation of the on-rung 
magnetization, where magnetization on odd rung is larger 
than on even rungs. This distribution remains almost un- 
changed within the plateau for < H < iJ + . 



3.3 Scaling properties of the magnetization plateau 



To find an accurate estimate on the critical exponent char- 
acterizing width of the magnetization plateau on the pa- 
rameter 6 we have computed the critical fields for the 
finite ladder systems with L = 6,8, 10, 12, 14, J± = 10 
and different values of the parameter S (0.01 < 5 < 0.1) 
and obtain their extrapolated values corresponding to the 
thermodynamic limit L — > oo. In Fig. [6] we have plotted 
the log-log plot of the plateau width versus 5. Calculations 
has been performed both in the case of ladder with antifer- 
romagnetic and ferromagnetic legs. We found that the best 
fit to our data (using the equation i7+ — ^ d'^) yields 
1/ = 1.82 ± 0.02 and v = 0.87 ± 0.01 for the ladders with 
ferromagnetic and antiferromagnetic legs respectively. 
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Fig. 6. Width of the magnetization plateau as a function of 
parameter i for 0.01 < <5 < 0.1 in the case of ladder with 
Jx = 10.0 and, (a) ferromagnetic (Jy = —1.0) and (b) antifer- 
romagnetic (J|| = 1.0) legs. 



4 Conclusion 

In this paper we have studied the elementary excitations 
and the magnetic ground state phase diagram of a spin 
5=1/2 two- leg ladder with alternating rung-exchange 
J±{n) = J±[l + (— 1)"5] using the continuum hmit bosoni- 
zation studies and the Lanczos method of numerical diag- 
onalizations for ladders up to L = 14. We have shown 
that the rung-exchange alternation leads to generation of 
a gap in the excitation spectrum of the system at mag- 
netization equal to the half of its saturation value. As 
a result of this new energy scale formation the magneti- 
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zation curve of the system M{H) exhibits a plateau at 
M = O.bMsat- The width of the plateau, is proportional 
to the excitation gap and scales as 5" , where critical expo- 
nent V = 0.87 ± 0.01 in the case of a ladder with isotropic 
antiferromagnetic legs and ly = 1.82 ± 0.01 in the case of 
a ladder with isotropic ferromagnetic legs. We have also 
calculated the magnetic field dependence of the on-rung 
spin-spin correlation functions. Comparison of these data 
for weak and strong rungs gave us an excellent description 
of the dynamics of the magnetization process in the case 
of ladder with non-equal rungs. 

In a standard way we estimate the magnetic condensa- 
tion energy as EmagiS)— Emag{0) ~ —(5^^. In the harmonic 
approximation the lattice deformation energy (per rung) 
can be estimated as Edef ~ S^. Therefore we can con- 
clude, that in the case of antiferromagnetic ladder, where 
2i/ = 1.64 < 2, the spontaneous appearance of an alter- 
nating rung exchange as a spin-Peierls instability during 
the magnetization process at M = O.bMsat is possible. 
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